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abstract 


This  report  describes  a  laser  Doppler  velocimeter  (LDV)  capable  of 
performing  measurements  in  vortex  cores,  presents  detailed  measurements 
°f  structure  and  stability  of  vortex  rings  and  also  presents  measure¬ 
ments  in  the  wake  of  a  model  helicopter  rotor. 

The  use  of  conventional  Instrumentation  such  as  hot  wire  anemometer 
for  measuring  velocity  in  vortex  filament  is  hampered  by  probe  inter¬ 
ference  and  the  difficulty  of  calibrating  the  probe  in  a  flow  field 
where  the  velocity  Is  changing  rapidly  in  magnitude  and  direction.  For 
these  reasons,  a  laser  Doppler  vel°cimeter  (LDV)  is  used  throughout 
these  experiments.  The  LDV  measures  the  velocity  of  smoke  particles  by 
measuring  the  Doppler  shift  of  laser  light  scattered  from  the  particles. 
The  LDV  system  designed  and  constructed  is  a  dual  scatter  system 
capable  of  determining  two  components  of  the  velocity  vector. 

The  most  unique  feature  of  the  LDV  system  is  the  electronics  system 
which  is  capable  of  handling  the  large,  fast  changes  in  velocity  which 
are  encountered  in  the  flow  fields  of  Interest.  The  electronics  system 
which  relies  on  the  fact  that  the  flow  field  is  periodic  requires  only 
a  specti  m  analyzer  for  demodulation  and  a  storage  oscilloscope  for 
retaining  the  data. 

A  detailed  study  of  the  structure  and  stability  of  vortex  rings  using 
the  two— component  LDV  system  has  been  made.  Vortex  rings  were  generate 
periodically  by  pulsing  air  through  a  sharp  edged  orifice  using  a  loud* 
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13.  ABSTRACT  (Continued) 

speaker.  The  circulation  and  core  size  a 
and  duration  of  the  electric  signal  appll 

Detailed  surveys  of  two  vortex  rings  we 
core  ring  with  a  ratio  of  core  radius  to 
relatively  thin  core  ring  with  a/r  -  .07 
that  the  vorticity  la  concentrated  in  a  s 
core  is  not  in  solid  body  rotation. 

During  the  course  of  the  investigation, 
are  unstable  was  rediscovered.  Azimuthal 
6  and  16  waves  spontaneously  develop  on  t! 
stable  and  the  smaller  the  vortex  core,  tj 
in  the  mode.  Theoretical  predictions  of 
but  the  correct  trend  is  indicated.  The 
ately  predicted  probahly  because  all  mode 
ficatlon  rate. 

A  single  bladed  model  helicopter  rotor 
blade  is  hollow  to  allow  smoke  to  be  depo 
desirable  point  along  the  span  providing 
terlng  particles  for  the  LDV  system. 

From  LDV  data  the  bound  circulation  and 
tribution  was  found.  Attempts  at  using  t 
structure  of  the  tip  vortex  were  largely 
fluctuations  In  the  position  of  the  vorte 
with  the  present  electronics  system. 
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CHAPTER  I 
INTRODUCTION 


1.1.  General  Comments 

The  current  interest  in  vortex  filaments  stems  from  the 
hazards  of  the  strong  vortex  wakes  of  jumbo  airliners  and  the 
noise  and  vibration  generated  by  blade  vortex  interactions  on 
helicopter  rotors.  Since  the  induced  motion  and  stability  of 
any  vortex  filament  depends  on  the  detailed  structure  of  the 
vortex  core,  it  is  desirable  to  make  experimental  measurements 
within  the  core.  The  purpose  of  this  thesis  is  to  describe  a 
laser  Doppler  velocimeter  (LDV)  capable  of  performing  measure¬ 
ments  in  vortex  cores,  to  present  detailed  measurements  of  the 
structure  and  stability  of  vortex  rings,  and  also  to  present 
measurements  of  the  flow  field  of  a  model  helicopter  rotor. 

1.2.  Literature  Review 

1.2.1.  Laser  Doppler  Velocimeter 

The  laser  Doppler  velocimeter,  a  device  that  measures  the 
velocity  of  particles  by  measuring  the  Doppler  shift  of  laser 
light  scattered  from  the  particles  was  introduced  in  1964  by 
Yeh  and  Cummins  (Ref.  1).  They  measured  the  velocity  of  poly¬ 
styrene  spheres  in  a  water  pipe  flow.  Since  then,  over  200 
papers  have  been  published  on  the  theory,  design  and  applica¬ 
tion  of  LDV  systems.  Two  bibliographies  of  LDV  literature  are 
available  (Refs.  2  and  3) ,  so  a  detailed  review  will  not  be 
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presented  here.  The  introduction  of  the  cross-beam  reference-beam 
system  by  Goldstein  and  Kreid  (Ref.  4)  in  1967  and  the  dual-scatter 
system  by  Rudd  (Ref.  5)  in  1969  have  brought  the  LDV  into  the 
realm  of  a  practical  fluid-measuring  device.  These  two  systems 
will  be  discussed  in  Chapter  II  where  the  relevant  equations  will 
be  derived  and  appropriate  references  cited. 

1.2.2.  VorLex  Rings 

The  study  of  vortex  rings  has  been  pursued  for  well  over 
a  century.  The  first  scientific  study  was  published  in  1858. 

W.  B.  Rogers  (Ref.  6)  observed  "rotating"  rings  in  water  and  air 
and  described  the  general  features :  the  generation  of  a  ring  by 
an  impulse  applied  across  an  orifice,  the  roll-up  of  the-  spiral 
sheets,  the  formation  of  a  translating  ring,  and  the  breakup  of 
rings  formed  from  droplets  of  a  heavier  fluid.  In  1867  W. 

Thomson  (Lord  Kelvin)  (Ref.  7)  in  a  note  to  1  ait's  translation 
of  Helmholtz's  paper  on  vortex  motion,  gave,  without  proof,  the 
correct  formulae  for  the  induced  velocity  of  a  ring.  He  then 
went  on  to  propose  (Ref.  8)  his  theory  of  vortex  atoms  based  on 
his  observations  of  the  indestructibility  of  vortex  rings. 

Reynolds  (Ref.  9)  in  1876  experimentally  observed  vortex  rings  and 
by  comparing  the  energy  and  momentum  of  the  ring  with  that  of  a 
circular  disk  he  derived  formulae  for  the  ring  and  core  radius 
relative  to  the  disk  radius.  In  his  monograph  (Ref.  10),  J.  J. 
Thompson  derives  formulae  for  the  induced  motion  and  modes  of  os¬ 
cillation  of  rings  which  are  valid  only  to  0(£n  a/rQ) .  These 
results  showed  the  rings  to  be  stable.  As  will  be  discussed  later, 
when  the  theory  is  carried  to  0(1)  the  vortex  ring  in  an  ideal 
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fluid  is  in  fact  unstable.  Hicks  (Ref.  11),  Dyson  (Ref.  12)  and 
Pocklington  (Ref.  13)  carried  the  theoretical  studies  of  the 
induced  velocity  to  higher  order  for  the  cases  of  solid  body 
rotation  cores  and  hollow  cores.  In  1894  Hill  (Ref.  14) 
derived  a  closed  form  solution  for  a  sphere  of  vorticity  whose 
strength  varies  linearly  with  distance  from  the  axis  of  symmetry. 
Northrop  (Ref.  15) ,  Sadron  (Ref.  16) ,and  Banerji  and  Barave 
(Ref.  17)  present  numerous  photographs  and  curves  of  distance 
traveled  versus  time,  for  rings  generated  by  different  impulses. 

In  a  rather  remarkable  paper,  Krutzsch  (Ref.  18)  experimentally 
observed  vortex  ring  instability.  His  results  showed  that  a 
sinusoidal  perturbation  (with  5  to  12  waves  around  the  perimeter) 
spontaneously  appears  on  the  filament  and  grows  in  amplitude. 
Krutzsch  also  measured  velocity  distribution  in  the  rings  using 
particle  tracer  techniques.  Although  Turner's  (Ref.  19)  main 
concern  was  buoyant  vortex  rings,  the  entrainment  ideas  he  pre¬ 
sents  are  also  relevant  to  non-buoyant  rings.  O'Brien  (20)  has 
derived  some  solutions  analogous  to  Hill's  spherical  vortex  for 
steady  spheroidal  vortices.  Magarvey  and  MacLatchy  (Ref.  21) 
have  some  nice  photographs  and  Zaroodny  (Ref.  22)  gives  a  review 
of  the  state  of  vortex  rings  up  to  1966.  Using  a  method  of 
matched  asymptotic  expansions,  Tung  and  Ting  (Ref.  23)  attempted 
to  obtai,'  the  induced  velocity  of  a  vortex  ring  with  a  viscously 
decaying  core  but  these  results  were  later  shown  to  be  in  error. 
Widnall,  Bliss,  and  Zalay  (Ref.  24),  Fraenkel  (Ref.  25),  and 
Saffman  (Ref.  26)  have  all  obtained  formulae  for  the  induced 
motion  of  small  core  rings  with  arbitrary  vorticity  distributions. 
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Johnson  (Ref.  27)  presents  a  summary  of  the  distance  vs.  time 
data  and  also  attempted  to  measure  the  structure  of  the  ring 
using  a  hot-wire  anemometer  but  due  to  interference  and  calibra¬ 
tion  problems  was  unsuccessful.  Recently  two  papers,  cited  in 
Refs.  28  and  29,  have  presented  experimental  flow  visualization 
results  for  rings  generated  by  distorted  holes  which  show 
various  instabilities.  Maxworthy  (Ref.  30)  also  observed  some 
instabilities  in  his  study  of  the  structure  of  rings.  However, 
his  attempts  to  measure  local  velocities  in  the  ring,  using 
hydrogen  bubble  techniques,  were  unsuccessful.  Finally,  Fraenkel 
'Ref.  31)  presents  some  examples  related  to  his  earlier  paper, 
including  a  higher  order  solution  for  a  ring  with  a  solid  body 
rotation  core. 

These  studies,  motivated  by  the  fundamental  interest  in  the 
motion  of  a  curved  vortex,  have  included  numerous  flow  visuali¬ 
zation  experiments  capable  of  measuring  the  ring  speed  and  ring 
size.  There  has  been  no  detailed  study  of  the  structure  of  the 
ring.  In  fact,  due  to  the  sensitivity  of  the  ring  to  probe 
interference,  even  the  total  circulation  of  the  ring  has  not  been 
measured.  In  the  present  study,  the  circulation,  stream  lines, 
and  vorticity  distribution  are  found  throughout  the  ring. 

1.2.3.  Helicopter  Rotors 

In  recent  years,  complex  free  wuhe  models  for  helicopter 
rotor  wakes  in  both  hover  and  forward  flight  have  been  developed 
by  Scully  (Ref.  32)  and  Langrebe  (Ref.  33).  Due  to  a  lack  of 
experimental  information,  various  assumptions  about  the  struc¬ 
ture  of  the  tip  vortex  and  its  subsequent  interaction  with  the 
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following  blade  were  made.  Ideally,  a  set  of  experiments  that 
measures  the  tip  vortex  structure  and  location,  and  detailed 
blade  loading,  is  needed.  The  structure  of  the  tip  vortex  for 
a  fixed  wing  in  a  wind  tunnel  has  been  measured  by  Chigier  and 
Corsiglia  (Ref.  34)  and  Rorke,  Moffitt  and  Ward  (Ref.  35) .  Their 
results  are  in  qualitative  agreement  but  differ  in  quantitative 
detail  indicating  the  need  for  more  extensive  investigation. 

Cook  (Ref.  36)  has  recently  made  some  preliminary  measurements 
in  the  tip  vortex  of  a  full-scale  rotor  on  a  test  stand.  His 
results  indicate  a  significantly  smaller  core  relative  to  the 
blade  chord  than  that  measured  by  Refs.  34  and  35. 
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CHAPTER  II 

THE  LASER  DOPPLER  VELOCIMETER 

2.1.  Introduction 

The  principle  of  operation  of  two  types  of  laser  Doppler 
velocimeters  will  be  considered  in  this  chapter.  A  detailed  des¬ 
cription  of  the  scattering  process,  operation,  construction,  and 
associated  electronics  for  one  of  these — the  dual  scatter  sys¬ 
tem —  will  also  be  given. 

2.2.  Basic  Principles  of  Laser  Doppler  Velocimetry 

Since  the  introduction  of  the  laser  Doppler  velocimeter 
(LDV)  in  1964  by  Yeh  and  Cummings  (Ref  1)  numerous  systems  have 
been  reported.  Two  of  these,  the  cross  beam- ref erence  beam  sys¬ 
tem  (Fig.  la)  of  ref.  2,  and  the  dual  scatter  system  (Fig.  lb) 
of  refs  3  and  37  will  be  considered.  These  two  were  chosen  because 
they  are  self  aligning;  that  is,  once  the  beams  are  inter¬ 

sect  at  the  probing  point  and  pass  through  the  aperture  in 
front  of  the  photomultiplier  tube,  the  system  is  aligned.  This 
eliminates  the  tedious  alignment  and  vibration  problems  of  other 
systems . 

Laser  light  that  is  scattered  from  particles  (smoke,  dirt, 
polystyrene  spheres)  moving  at  the  local  fluid  velocity  under¬ 
goes  a  Doppler  shift  given  by 

+0(i|9  U> 
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where 

A 

f  ,  X,  kj  *  frequency,  wavelength  and  unit  vector  in  the 
propagation  direction  of  the  incident  laser 
beam 

V  =  velocity  of  particle 

A 

k  =  unit  vector  in  the  direction  that  the  scattered 

light  is  detected 

c  =  velocity  of  light 

The  frequency  of  the  scattered  radiation  f,  which  is  of  the 

14 

order  of  light  frequencies  (^5  x  10  cy/sec)  can  be  detected 
either  directly  using  a  Fabry-Perot  interferometer  or  can  be 
heterodyned  with  some  light  wave  at  a  known  frequency  to  obtain 
a  signal  whose  frequency  is  equal  to  the  shift  of  the  laser  light 
frequency.  The  present  resolution  of  Fabry-Perot  interferometers 
is  only  ^2MHz  which  implies  a  velocity  resolution  of  ^10  ft/sec. 
For  the  present  study  this  is  inadequate  so  a  heterodyning  system 
is  used. 

In  the  two  types  of  heterodyne  systems  studied  the  cross¬ 
beam  reference  beam  system  (Fig.  la)  heterodynes  the  scattered 
radiation  with  radiation  at  the  laser  frequency  fQ  and  the 
dual  scatter  system  (Fig.  lb)  heterodynes  the  scattered  radiation 
of  two  different  incident  beams.  In  both  cases  the  difference 
frequency  is  calculated,  using  Eq.  1  and  the  geometry  of  Fig.  1, 
to  be 


2.  U  SfN 

^ C* 
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where  U  is  the  velocity  component  in  a  direction  perpendicular  to 
the  bisector  of  the  angle  between  two  beams.  Fig.  2  shows  fD  ver- 

O 

sus  U  for  various  values  of  0  with  A=6328A,  the  wavelength  of  a  He 
Ne  laser.  It  should  be  noted  that  there  is  a  180°  amgiguity  in 
determining  direction  of  U.  The  plus  or  minus  direction  of  U  must 
be  determined  by  prior  knowledge  or  flow  visualization. 

The  output  signal  of  the  photomultiplier  (P.M.)  tube  is  pro¬ 
portional  to  the  intensity  which  for  the  reference  beam  system  is 
***  =  (£«♦  **)■(£*+ 

=  i-ZEK-E%(t)cos(ZTTf9t) 

Er  =  electric  field  of  reference  beam 
^  =  scattered  electric  field 


when  a  single  particle  traverses  the  measuring  volume.  This  pro¬ 
duces  an  output  like  that  shown  in  Fig.  3a.  The  d.c.  level  is 

essentially  constant  since  E  (t)<<E  .  The  a.c.  component  is  modu- 

s  c 

lated  by  the  term,  EE  (t)  which  is  generally  a  gaussian. 

K  S 

In  the  dual  scatter  system  the  light  scattered  from  one  beam  is 
heterodyned  with  light  scattered  from  the  other  beam  resulting  in 

x„,-  i- 1  *  2  cos(j-rr  v) 

=  electric  field  scattered  by  beam  one 


£5  =  electric  field  scattered  by  beam  two 

The  results  for  E  (t)  =  E  (t)  are  shown  in  Fig.  3b.  The  d.c. 

S1  s2 

component  is  now  replaced  by  a  low  frequency  component.  The  actual 
signal  when  a  .0002  wire  traverses  the  measuring  volume  is  shown  in 
the  oscilloscope  trace  of  Fig.  4a.  When  the  wire  doesn't  pass 


5 


i 
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through  the  center  of  the  measuring  volume  a  signal  (Fig.  4b)  con¬ 
taining  a  large  amount  of  the  low  frequency  component  and  small 
amount  of  the  Doppler  signal  component  is  obtained. 

In  the  present  investigation  more  than  one  scatterer  will 
be  in  the  measuring  volume  at  one  time  resulting  in  a  signal  which 
is  the  sum  of  the  signals  from  individual  particles.  A  typical 
"continuous"  signal  for  a  dual  scatter  system  is  shown  in  Fig.  4c. 
This  signal  was  obtained  from  a  rotating  ground  glass  disk.  It 
shows  the  general  features  of  dual  scattering;  namely,  a  low  fre¬ 
quency  component  and  low  frequency  modulation  of  the  Doppler 
signal.  The  low  frequency  modulation  is  a  result  of  random 
constructive  and  destructive  interference  among  the  waves  scat¬ 
tered  by  many  particles,  causing  the  signal  to  go  randomly  to 
zero  during  severe  destructive  interference.  This  phenomenon, 
called  "drop-out",  and  the  low  frequency  component,  are  the  major 
problems  in  the  electronic  processing  of  the  signal. 

A  comparison  of  the  frequency  spectrums  of  the  two  systems 
is  shown  in  Fig.  5  for  scattering  from  a  ground  glass  disk  and  a 
turbulent  jet.  Although  the  magnitudes  of  the  signals  are  about 
equal  the  low  frequency  components  in  the  dual  scatter  system 
are  greater. 

The  object  of  this  investigation  is  to  determine  which  of 
the  two  systems  will  give  a  better  signal-to-noise  ratio.  The 
reference  beam  system  has  the  advantage  that  the  signal  can  be 
increased  by  increasing  the  strength  of  the  reference  beam. 
However,  the  solid  angle  through  which  scattered  light  may  be 
collected  is  limited  since  the  Doppler  shift 


V.-  ( h,  -  hi) 
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is  a  function  of  viewing  direction  kg  so  that  at  large  solid 
angles  the  Doppler  frequency  would  be  unacceptably  broadened. 

.  A  -A  % 

On  the  other  hand,  the  frequency  ~*  (  ^1*  of  a  dual 

A 

scatter  system  depends  only  on  the  incident  beam  directions 

A  A 

k  ,k  so  that  the  solid  angle  for  viewing  scattered  radiation 
1  2 

is  unlimited.  However,  the  strength  of  the  signal  from  a  dual 
scatter  system  is  determined  entirely  by  the  scattering  process. 
Also  there  will  be  more  low  frequency  components  in  the  dual 
scatter  system. 

This  initial  investigation  did  not  indicate  that  one  sys¬ 
tem  was  clearly  superior  but  that  both  performed  adequately. 

At  this  point  it  was  felt  that  further  comparison  would  not  be 
fruitful,  and  the  dual  scatter  system  was  chosen  for  further 
study. 

2.3.  Scattering  By  a  Spherical  Particle 

The  far  field  solution  for  the  electric  field  of  an  inci¬ 
dent  plane  wave  traveling  in  the  z-direction  and  linearly 
polarized  in  the  x-direction 

-►  ^  r  t  (u;-t  -  kSL) 

£  =  =  incident  electric  field 


scattered  by  a  dielectric  sphere,  of  diameter  d,  located  at  the 
origin  of  the  coordinate  system  (Fig.  6)  is  given  by  the  Mie 
solution  (Ref.  38).  N 

E  =•  i  f-sm  <^>  S^C©,0^)  e* 

s  k  r  u 


+  cos 


4>  Sl,C©.0tr>)  £«.] 
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The  complex  scattering  amplitude  Sxis  the  scattered  component 
whose  electric  vector  is  perpendicular  to  the  plane  of  observa¬ 
tion  (the  plane  determined  by  the  direction  of  observation  and 
the  direction  of  propagation  of  the  incident  beam)  and  the 

component  parallel  to  the  plane  of  observation.  The  amplitudes 
Sx  and  S„  are  functions  of  the  azimuthal  angle  0,  the  size 
parameter  a 

oc-  -  TTfll 

r-  ~  X 

and  the  index  of  refraction  of  the  sphere,  n.  The  solution 

is  a  spherical  wave  whose  polarization  depends  on  the  relative  amp- 

A  A 

litudes  and  phases  of  the  e^  and  e0  components. 

The  intensity  of  the  scattered  radiation  is 

X  =  +  i„cosr<p] 

where  LA  =  ISJ*  <- „  =  I  ^>..1^ 

The  real  and  imaginary  parts  of  Sx  and  are  shown  in  Fig.  7 

•  • 

and  tx  and  L„  are  shown  in  Fig.  8.  Note  that  for  small  0 

•  • 

the  values  of  lx  and  l,,  are  approximately  equal  and  also  that 

hhorp  i  q  nr*  nhaca  chi  ff  hofuoon  hhom  f  i  o  ^  W  ( ^  ( ^l|\  ) 


there  is  no  phase  shift  between  them  (i.e., - v  *  ■  —  ] 

This  implies  that  the  scattered  radiation  is  linearly  polarized. 
For  small  values  of  0  the  complex  amplitudes  Sx  and  Sj| 


are  equal  so 


P  ^  C 

*  kr 


c(wt-fcr) 


<5^-  sm$  cos  <$> 
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If  a  lens  of  focal  length  f  is  aligned  on  the  z-axis  the  scattered 
radiation  is  transformed  to  e  plane  wave 


.  C(u»t-k€) 

L  Co 


e 


* 


The  polarization  of  this  wave  is  the  same  as  the  incident  polari¬ 
zation,  a  fact  that  will  be  used  later  for  separating  the  com¬ 
ponents  of  the  velocity  vector.  Note  that  this  is  only  valid 
for  small  0,  and  that  the  polarization  for  large  0  will  not  in 
general  be  the  same  as  the  incident  wave. 

In  the  actual  LDV  system,  the  laser  beam  is  focused  by  a 
lens  of  focal  length  f.  The  electric  field  in  the  focal  region 
is  given  by  the  plane  wave  (Ref  46) 


C-x  G 


W)  x  radius  of  focused  spot  to  —  intensity  points 

•  IT  e2 

_  Sf.  Focal  length 

1  Diameter  of  unfocused  laser  beam(l/e2  intensity) 

and  I0  is  the  total  integrated  intensity.  The  power  of  the  laser 
is  related  to  I0  by 


I0  = 


2P. 


TTb  ' 


which  was  obtained  by  integrating  E  •  E  over  the  x-y  plane. 

Although  the  magnitude  of  the  electric  field  varies  over 
the  particle,  this  effect  will  be  small  since  the  laser  spot  size 
will  be  ^100  times  the  diameter  of  the  particle.  The  electric 
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field  scattered  by  a  particle  in  the  vicinity  of  a  focused  laser 
will  be 

-  ;  Ct  VA  '  t <.*■>•«: -hr) 

4  Tr  e  x 


For  future  use  the  scattered  field  from  a  particle  located 

-f  4 

at  an  arbitrary  point  r^  by  a  plane  wave  with  vector  k^  (lying 
in  z,y  plane)  is  found  to  be 

^<4  e;  Co**  -  K^l  -  T;  •  rr] 


Kinl 


£  -  ■sin  4>i.  Sj.Ce^.n)  e  +  c.os  <j>i  £»r} 


where  the  subscripted  coordinates  refer  to  the  coordinate  system 
shown  in  Fig.  9.  Since  | R | > > | r  j  the  origin  of  this  coordinate 
system  (x^y^z^  can  be  taken  at  the  origin  of(x,y,z)  so  that  z^ 

A 

lies  along  k^  and  y^  is  in  the  z,y  plane.  Also 

h|nl-  k(K-rP(s  h(|R|-|--v)£  **'*»'**• 
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So  the  scattered  field  is  approximately 

~i\C*P£  * 

F.  ^  >  (x«.)  <=•  e.* 

*  HR 

[,-%tw4>i  S^CQj ,«K,n)e4.  +  co%^£S1,(oi«x.n)e.r] 

(2) 

2.4.  Detailed  Description  of  Dual  Scatter  System 

In  the  dual  scatter  system,  two  parallel  laser  beams  are 
focused  at  a  common  measuring  point  as  shown  in  Fig.  10.  The  ra¬ 
diation  scattered  by  beam  No.  1  is  given  by  Eq.  2  with  i  ■  1  and 
that  scattered  by  beam  No.  2  is  also  given  by  Eq.  2  with  i  *  2. 

The  subscripted  coordinates  (x^y^z^)  or(r^,  6^,  refer  to 
coordinates  whose  z^  axis  lies  along  with  both  Jc^  and  y^  in 
the  z,y  plane.  Since  the  scattering  is  linear  the  electric  field 
at  any  distant  point  is  just  the  sum  of  the  two  scattered  waves. 

The  intensity  of  light,  I,  scattered  by  the  two  beams  at  any 
point  on  a  sphere  of  radius  R  is  given  by 


f  -  £?(.*?*  y\  co»l  tin1  •»') 

'rSv  e  * 

L  o',  +  a-t  t  to;,  cos(  O) 


where 


=  I-.  e 


L  •»>  im1  <#>,  LJ_(e.')  +  t.|i 


-  I»_  e 


+  ^  L  «w  &J 


[  SIM*  4>x  ^(©t)  +  tOS*4>t  i.,  C©t^] 


(r<.,I<0  4>,  s iw  <k  (tj.Ce;)  tx (®J)  e+i-  e^t 

+  cos  <f>,  L°^  ^>t  (  t „ (©^  t„ ( e*))^  €^  e-ot 

+  sim  <p,  cos  <£t  ( txc©^  i„ce^)4  e.,-e*t 

+  cos  <#>,  SIM  <j>t(  tj.(©^  e.t-e^3 


Several  conclusions  can  be  drawn  from  Eq.  (3)  .  Setting  ^"p  “  ^op  +* 

C  *<>,.♦  «t ,  *  vt ,  z„f  +  wt ) 

(valid  to  Jli)  the  a.c.  component  has  a  frequency 

c  z  siw  Ss  lor! 
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i.e.  ,  by  measuring  this  frequency  and  0g  the  component  of  velocity 
in  the  plane  of  the  two  beams  and  in  a  direction  perpendicular 
to  the  bisector  of  two  beams  can  be  determined.  The  exponent 
indicates  that  when  a  particle  is  on  the  ellipsoid 


V  b./co*e%  ) 


the  intensity  will  be  1/e^  of  the  intensity  for  a  particle  located 

at  (0,0,0).  This  ellipsoid  is  called  the  probe  volume.  Although 

the  probe  volume  is  actually  determined  by  the  minimum  detectable 
2 

signal  the  1/e  points  provide  a  good  estimate.  The  volume  of 


the  probe  is 


3  sin8. 


For  small  values  of  0B ,  the  probe  is  approximately  a  cylinder 
with  radius  bQ  and  length  2bQ/sin  (0g/2). 

The  signal  from  a  single  particle  given  by  Eq.  3  is  the 
sum  of  two  gauss ian  pulses  and  the  Doppler  shift  term  with  the 
resulting  signal  discussed  earlier  and  shown  in  Fig.  4a.  The 
length  of  tne  Doppler  pulse  is  controlled  by  a  gaussian  envelope 

cos5®}  +(io, 


The  wider  this  gaussian  envelope  the  less  broad  will  be  the  sig¬ 
nal  in  the  frequency  domain.  The  best  case  occurs  for  V  =  W  =  0 

x  =  y  =  z  =0  giving  a  width 
°f  °P  °P 

Af  1  0R/, 

T  ~  2f  cot  B'2 
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0 


It  should  be  noted  that  this  is  the  minimum  broadening 
that  can  be  expected.  Any  particle  with  a  V  and/or  W  component 
and/or  not  passing  directly  through  the  center  of  the  scatter¬ 
ing  volume  (x  /0,  y  ,  zQ  /0)  will  cause  additional  broad- 
P  P  P 

ening. 

A  useful  approximation,  and  one  that  is  often  valid  in 
practice,  is  to  assume  8fi<<l  (or  (it-  0B)<<1)  and  the  solid  angle 
(ft)  over  which  scattered  radiation  is  collected  is  also  small. 


Then 


£  4>*  *  <t> 

e,  a  e 

i-x  ( e,)  -  tx(e) 


Also  letting  IQ  =  IQ  =IQ  which  is  the  optimum  case. 


X  -  2  £a(®)  *  co 


[coSH(y^r^)  +  COS2Tr(^3U"' 


The  term  in  braces  is  simply  the  intensity  distribution  in  the 
probe  volume.  This  intensity  distribution  oscillates  from  4Iq 
to  0  thus  giving  rise  to  a  set  of  interference  fringes.  The 
fringe  model,  first  suggested  by  Rudd,  implies  that  as  a  par¬ 
ticle  moves  into  a  bright  fringe,  it  scatters  light,  then  it 
moves  into  a  dark  fringe,  then  back  into  a  bright  fringe 


scattering  light,  etc.  Note  that  this  model  is  only  valid  for 
small  forward  or  backward  scattering. 

Until  this  point,  the  polarization  of  the  incident  beams 
has  been  in  the  x-direction,  i.e.,  perpendicular  to  the  plane 
of  the  two  beams.  For  the  small  angle  scattering  case  this 
maximizes  the  Doppler  shift  term  of  the  intensity  since  the 
fringe  contrast  will  be  maximum.  To  show  this  consider  the  a.c. 
part  of  the  intensity  within  the  probe  volume  formed  by  two 
incident  beams  with  arbitrary  linear  polarization 

Ex,  =  Eo.,  Eb, 

A  A 

where  the  unit  vectors  a,  b  are  shown  in  Section  AA  of  Fig.  1. 
The  a.c.  or  Doppler  shift  terra  is 

c.o?,e6 +Eb>Et>J  (4> 

For  arbitrary  polarization  this  equation  shows  that  the  signal 
decreases  as  0B  increases.  If  incident  beams  are  polarized  in 

A  /\ 

the  b.  and  b~  directions  (i.e.  E,  =  E  =  =  0  )  then  the  maximum 

J.  2  ax  a2 

signal  will  be  obtained.  In  cases  where  the  scattering  is  not 
limited  to  small  angles  the  polarization  of  the  scattered  waves 
(i.e.  Sx  and  )  must  be  considered. 

2.5.  Two  Dimensional  LDV  System  and  Associated  Electronics 

The  LDV  system  constructed  for  extensive  measurements  is 


a  dual  scatter  system.  A  one-dimensional  dual  scatter  system  is 
shown  in  Fig.  11.  The  laser  beam  is  divided  by  a  standard  prism- 
type  beamsplitter  with  a  non-absorbing  beamsplitting  interface. 
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This  method  of  dividing  the  laser  beam  is  preferable  to  other  sys¬ 
tems  using  mirrors  and  beamsplitters  since  only  one  standard 
optical  component  is  needed  and  the  path  lengths  (distance  from 
laser  to  measuring  point)  of  two  beams  are  equal. 

A  two-dimensional  system  is  constructed  by  adding  a  second 
beamsplitter  (rotated  90  degrees  with  respect  to  the  first)  to 
form  four  parallel  beams.  The  pattern  of  the  four  beams  is  shown 
in  Fig.  12.  The  polarization  of  the  beams  is  chosen  to  maximize 
the  signal  (Eq.  4)  and  to  eliminate  cross  talk  between  the  two 
dimensions.  By  placing  an  analyzing  polaroid  in  front  of  the  P.M. 
tube  either  of  the  velocity  components  can  be  chosen  with  no 
interference  from  the  other  component.  This  is  possible  because 
the  polarization  is  changed  only  slightly  by  the  scattering  pro¬ 
cess.  In  the  present  system,  both  components  are  processed 
simultaneously  through  one  P.M.  tube  using  the  electronic  system 
described  below.  A  photograph  of  the  two-dimensional  LDV  system 
is  shown  in  Fig.  13.  The  two  components  could  be  processed 
separately  by  splitting  the  scattered  radiation  with  a  polarizing 
type  beamsplitter  and  using  two  P.M.  tubes  and  two  sets  of  elec¬ 
tronics. 

The  electronic  system  for  processing  the  LDV  signal  is  one 
suggested  by  Ref.  39  for  use  with  periodic  flow  fields.  This 
"sampling"  system  will  be  described  for  use  with  periodically 
produced  vortex  rings  in  terms  of  the  time  diagrams  of  Fig.  14 
and  the  equipment  layout  of  Fig.  11.  The  function  generator 
initiates  a  pulse  which  is  amplified  and  used  to  drive  a  loud¬ 
speaker  which  forms  a  vortex  ring  at  a  sharp  edged  orifice.  P 
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synchronous  pulse  is  also  taken  off  of  the  function  generator 
(Fig.  14a) . 

The  velocity  versus  time,  as  the  vortex  ring  goes  past  the 
measuring  point,  is  shown  in  Fig.  14b.  This  is  proportional  to 
the  frequency  versus  time  of  the  output  signal  of  the  P.M.  tube 
which  is  fed  into  the  input  of  the  spectrum  analyzer.  The  syn¬ 
chronous  pulse  from  the  function  generator  initiates  a  time  delay 
which  starts  the  sawtooth  sweep  generator  (Fig.  14c)  for  the 
spectrum  analyzer  after  a  time  tp.  The  linear  sawtooth  repre¬ 
sents  the  center  frequency  of  the  narrow  filter  of  the  spectrum 
analyzer.  When  the  frequency  from  the  P.M.  tube  is  equal  to  the 
instantaneous  frequency  of  the  narrow  band  pass  filter  of  the  spec¬ 
trum  analyzer  a  "pip"  will  occur  on  the  spectrum  analyzer  face 
and  at  the  spectrum  analyzer  output  (Fig.  14d) .  This  condition 
can  be  determined  by  superposing  Fig. 14c  on  Fig.  14b.  The  linear 
sawtooth  is  fed  to  the  input  of  a  gate  which  is  only  opened  when 
a  "pip"  occurs  on  the  spectrum  analyzer.  The  output  of  the  gate 
—  in  this  case,  three  spots  —  is  stored  on  a  storage  oscillo¬ 
scope.  The  three  spots  correspond  to  a  zero  frequency  marker  and 
two  points  on  the  velocity  versus  time  curve  we  wish  to  determine. 
By  varying  the  time  delay  the  remaining  points  on  the  velocity 
curve  should  be  stored. 

The  results  of  demodulating  a  test  oscillator  signal  varying 
from  50  kHz  to  950  kHz  at  a  rate  of  33  Hz  are  shown  in  Fig.  15. 

The  entire  LDV  system  was  checked  by  measuring  two  components 
of  the  velocity  of  a  spinning  disk  as  the  disk  was  traversed 
across  the  measuring  point.  The  x-component  of  velocity  U  =  -  toy 


36 


should  be  a  constant,  and  the  y-component  V  =  ojx  should  vary  lin¬ 
early  as  the  disk  moves  in  the  x-direction.  The  results  of  this 
check,  shown  in  Fig.  16,  indicate  that  the  linearity  of  the 
electronics  system  is  excellent  and  that  the  velocity  error  is 
within  the  error  in  determining  the  r.p.m.  of  the  disk.  The 
ambiguity  in  determining  the  direction  of  U  and  V  is  demon¬ 
strated  in  Fig.  16.  The  U-comporient  is  always  negative  and  the 
V-component  is  negative  for  negative  values  of  X. 

The  two-dimensional  system  has  been  used  to  measure  the  velo¬ 
city  distribution  in  periodically  produced  vortex  rings  (smoke 
rings) .  (The  results  of  extensive  measurements  are  presented 
in  Chapter  IV.)  The  rings  are  similar  to  a  helicopter  rotor 

wake  since  regions  of  time  varying  concentrated  vorticity  are 
produced,  yet  the  flow  field  is  simpler  since  it  is  axisymmetric. 

A  picture  of  a  vortex  ring  is  shown  in  Fig.  17  and  the  correspond¬ 
ing  velocity  distributions  in  Fig.  18.  The  velocity  curves  are 
the  U  and  V  components  of  velocity  versus  time  as  the  vortex  ring 
travels  by  the  fixed  measuring  point.  There  is  a  180°  ambiguity 
in  the  velocity  vector  (since  frequency  is  always  positive)  .  Flow 
visualization  indicates  that  the  U-component  is  always  positive 
while  the  V-component  is  positive  to  the  left  of  the  central  zero 
and  negative  to  the  right. 

2.6.  Summary  of  the  LDV  System 

Using  the  formulae  derived  above,  the  characteristics  of  the 
LDV  system  for  a  half  angle  of  0g/2  =  3.90°  are  calculated  and 
summarized  in  Table  1.  For  this  angle,  which  is  typical  of  all 


V 


37 

systems  used,  the  volume  from  which  adequate  Doppler  signal  is 
received  is  aporoximately  a  cylinder  .0026"  in  diameter  and 
.038"  long.  On  the  average  there  are  about  five  l.Oy  diameter 
particles  in  this  cy1inder  at  one  time. 
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CHAPTER  III 

PARTICLES  FOR  USE  WITH  LDV  SYSTEM 

3.1.  Introduction 

A  laser  Doppler  velocimeter  measures  the  velocity  of  parti¬ 
cles  which  are  generally  assumed  to  have  a  velocity  equal  to  the 
local  fluid  velocity.  Although  it  is  theoretically  possible  to 
calculate  the  fluid  velocity  from  the  particle  velocity,  if  they 
are  not  equal  this  would  require  detailed  knowledge  of  instan¬ 
taneous  particle  size  which  is  usually  not  available.  Thus  it 
is  important,  especially  for  measurements  in  air  where  the  par¬ 
ticle  density  is  several  orders  of  magnitude  times  the  fluid 
density,  to  assure  that  the  particles  are  small  enough  to  follow 
the  flow. 

The  purpose  of  this  section  is  to  discuss  a  method  of  gen¬ 
erating  suitable  particles,  and  measuring  their  diameter  and  num¬ 
ber  density.  Also  a  procedure  for  estimating  the  difference 
between  the  fluid  and  particle  velocity  will  be  described. 

3.2.  Particle  Generation  and  Measurement 

In  the  present  experiment  oil  particles  were  generated  by 

* 

a  Farval  vortex  mist  lubricator  ,  a  device  normally  used  for 
lubricating  high  speed  bearings.  The  oil  particles  are  generated 
by  mechanically  breaking  up  oil  drops  at  a  sonic  orifice.  This 

* 

Farval  Mist  Lubricator  Model  No.  VM30B .  Available  through  Fluid 
Power  Division  of  Eaton,  Yale,  and  Towne,  Inc.  of  Cleveland,  Ohio. 
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results,  after  filtering,  in  cold  oil  particles  in  the  micron 
diameter  size  range.  This  method  of  obtaining  particles  was  found 
to  be  extremely  convenient,  since  it  is  only  necessary  to  turn 
on  the  'hop  air  to  obtain  suitable  particles. 

* 

The  oil  used  in  the  lubricator  is  called  "Fog  Juice"  ,  a 
light  oil  used  in  theatrics  for  producing  smoke  for  graveyard 
scenes,  etc.  This  oil  has  the  advantage  that  it  will  evaporate 
in  about  a  week,  so  the  laboratory  and  equipment  do  not  become 
covered  with  oil  film. 

The  lubricator  manufacturer  specifies  that  the  largest 
particle  size  is  approximately  2\i  in  diameter  but  the  size  is  a 
function  of  the  oil  used.  In  order  to  determine  the  size  more 
accurately  the  diameter  was  found  by  measuring  the  intensity  dis¬ 
tribution  in  the  forward  scattering  lobe.  This  technique, 
developed  by  Hodkinson,  Ref  40,  uses  the  fact  that  the  relative 
angular  distribution  of  scattered  intensity  within  the  Mie 
forward  lobe  is  a  function  of  size  only  and  not  refractive  index. 

The  scattering  function  was  measured  from  5°  to 

150°  using  the  setup  shown  in  Fig.  19.  By  forming  the  ratio  of 
tjL  at  various  angles  the  value  of  a  =  ^  can  be  found  using 
Fig.  20,  which  is  taken  from  Hodkinson. 

The  scatter  in  the  data  shown  in  Fig.  20  indicates  that 
there  is  a  distribution  of  particle  sizes  but  the  mean  size  is 
a  =  5.0  which  gives  a  mean  diameter  of  d  =  l.Op. 


*  Type  1964  Fog  Juice,  available  from  Mole-Richardson  Co., 
Hollywood,  Calif. 
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The  loss  in  intensity  of  a  beam  of  light  traversing  a  dis- 

3 

tance  A  through  a  dispersion  of  N  particles/cm  is 

x;=  e 

where  Q  is  the  scattering  efficiency  which  is  equal  to  th 

SCd 

2 

scattering  cross-section  divided  by  tt  ( d/2 )  .  By  placing  the 

detector  of  Fig.  II  at  0  =  0  and  measuring  the  ratio  of  inten¬ 
sities  with  and  without  particles  in  the  cell  the  number  of  par- 
ticles/cm'*  can  be  determined.  The  results  of  this  measurement, 
using  a  value  Q  =3.0  from  Ref  38  for  a  =  5.0  and  n  =  1.5, 

SCd 

6  3 

give  N  =  8.5  x  10  particles/cm  . 

3.3.  Motion  of  Spherical  Particles  in  a  Flow  Field 

In  order  to  estimate  the  difference  between  fluid  and 
particle  velocities  the  equation  of  motion  of  the  particle  is 
non-dimensionalized  and  the  appropriate  dimensionless  parameters 
derived. 

The  equation  of  motion  of  a  small  spherical  particle  sub¬ 
ject  to  Stokes  drag  is  (Ref  41) 


TT 

<o 


-  3TrV+ff^(Up-  «*) 
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where  U  =  velocity 
p  =  density 
v  =  viscosity 
d  =  diameter  of  particle 

Subscripts  p  and  f  refer  to  particle  and  fluid  respectively. 
The  meaning  of  the  various  terms  is  as  follows. 

The  term  on  the  left-hand  side  is  the  force  required  to 
accelerate  the  particle.  The  first  term  on  the  right-hand  side 
is  the  viscous  resistance  force  according  to  Stokes'  law.  The 
second  term  is  due  to  the  pressure  gradient  in  the  fluid  surround¬ 
ing  the  particle,  caused  by  acceleration  of  the  fluid.  The  third 
term  is  the  force  to  accelerate  the  apparent  mass  of  the  particle 
relative  to  the  ambient  fluid.  The  fourth  term  is  the  so-called 
"Basset"  term,  ich  takes  account  of  the  effect  of  the  deviation 
in  the  Stokes  flc*/  pattern  from  steady  state.  The  last  term  is 
the  force  due  to  gravity. 


Non-dimensionalizing  with  respect  to  a  reference  fluid  velocity 

UREF  and  len<?th  Lr5f  and  using  V  ^  =  "fa 

the  equation  of  motion  becomes 

-E>^)  +  ( £  -«o  -  C  f  zM') 

i-  G  ~  O 

where  all  quantities  are  nondimensional 

~  (  I  +  ^  =  Stokes  number 

't  =  ^  h. 

/s^  />* 
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In  the  present  study  oil  particles  are  used  in  air  so 


Pf  -3 

—  =10  implying  B  <<  1  and  that  C  <<  Stk  so  the  equation  can 

PP 


be  written  approximately  (The  gravity  force  is  also  neglected) 


k +(u/~  «*)  =  o 

(\t 


The  role  of  the  time  scale  t  can  be  seen  by  considering 
the  simple  problem  of  a  particle  moving  in  a  straight  line  with 


a  velocity  at  t  =  0  of  UQ  in  a  stationary  fluid.  At  any  time  t 


later  the  velocity  will  be 

CU 


u, 


^  e 


showing  that  t  is  the  characteristic  rise  time  for  a  particle 
to  reach  the  local  fluid  velocity. 

For  one  micron  diameter  oil  particle  in  air  r  =  4  x  10  6sec 
and  the  Stokes  number,  which  is  a  measure  of  the  difference  between 


fluid  and  particle  velocities,  is 
-3 


st; 

Thus 


c  =  4  x  10 


U 


REF 


for  =  10  ft/sec  and  L^p  =  .01  ft. 


4% . 
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CHAPTER  IV 
VORTEX  RINGS 

4.1.  Introduction 

Although  the  fundamental  interest  in  vortex  rings  has  led 
to  many  flow  visualization  studies,  there  have  been  relatively  few 
attempts  to  measure  the  detailed  structure  of  the  ring.  This  is 
due  in  part  to  the  sensitivity  of  flows  with  concentrated  vorti- 
city  to  the  insertion  of  a  probe. 

Tracer  techniques  were  attempted  by  Krutzsch  (Ref  18)  in 
1939  using  aluminum  flakes  and  Maxt/orthy  (Ref.  30)  in  1970  using 
hydrogen  bubbles  but  due  to  strong  axial  gradients  and  associated 
radial  velocities  their  results  were  inaccurate.  Hot  wire  measure¬ 
ments  (Ref  27)  are  hampered  by  probe  interference  and  the  diffi¬ 
culty  of  calibrating  the  probe  in  a  flow  field  where  velocity  is 
changing  rapidly  in  magnitude  and  direction.  These  problems  are 
circumvented  in  the  present  experiment  by  use  of  a  laser  Doppler 
velocimeter  (LDV) . 

Although  tne  formation  process  cannot  be  predicted,  there  are 
two  analytical  treatments  of  steady  vortex  rings,  the  small  core 
vortex  where  the  ratio  of  core  size  to  ring  radius  is  small  and 
Hill's  spherical  vortex  where  vorticity  is  distributed  throughout  a 
sphere.  Several  recent  studies  (Refs  24&26)  of  the  motion  of 
small  core  vortex  rings  give  the  translational  velocity  of  the 
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ring  as  (from  Ref  24) 

«o-  -yQ 

where  A  is  an  integral  over  the  tangential  velocity  distribution 
in  the  core  (defined  later) .  For  solid  body  rotation  in  the  core 
A  =  1/4  giving  the  classical  Kelvin  results.  From  the  above 
equation  the  relevant  non-dimensional  parameter  to  characterize 
the  structure  of  the  vortex  core  is 

V  *  r0 

Since  the  term  in  braces  depends  on  the  detailed  distribution  of 

vorticity  it  is  of  interest  to  measure  this  distribution.  Also, 

models  for  studying  the  stability  or  entrainment  will  depend  on 

the  detailed  structure  of  the  ring.  In  the  present  study  the 

circulation,  streamlinas  and  vorticity  distribution  are  found 

throughout  the  ring  from  measurements  of  the  axial  and  radial 

velocities  made  prior  to  instability  and  compared  to  theoretical 
predictions. 

4.2.  Description  of  the  Experiment 

Vortex  rings  are  generated  by  pulsing  air  through  a  sharp- 
edged  orifice  using  a  loudspeaker.  (The  circulation  and  core 
size  are  controlled  by  the  amplitude  and  duration  of  the  elec¬ 
trical  signal  applied  to  the  loudspeaker.)  (See  Fig.  11.) 

The  electronic  system  for  processing  the  LDV  signal  is  a 
data  sampling  system  that  requires  a  periodic  flow  field.  To 
obtain  a  complete  velocity  profile,  the  rings  are  produced 
periodically;  the  distance  between  rings  was  approximately  100 
ring  radii  to  prevent  any  interaction  between  the  rings. 
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The  LDV  system  measures  two  components  of  velocity  versus 
time  at  a  fixed  point  in  space  as  the  vortex  ring  goes  by  the 
measuring  point.  The  velocity  distribution  is  transformed  to 
the  steady  coordinate  system  of  Fig.  21  by  letting  Z  =  U  t  where 
UQ  is  the  instantaneous  velocity  of  the  vortex  ring  at  the  mea¬ 
suring  point.  This  quasi-steady  approximation  assumes  that  the 
properties  of  the  vortex  ring  do  not  change  (due  to  entrainment 
of  surrounding  fluid  or  wake  shedding)  signif icantly  in  the  time 
it  takes  the  vortex  ring  to  pass  the  measuring  point.  This  is 
justified  by  flow  visualization  studies  of  the  variation  of  ring 
speed  UQ  and  size  of  the  ring. 

Detailed  surveys  of  a  relatively  thick  core  ring  (denoted 
Ring  No.  1)  and  a  thin  core  ring  (Ring  No.  2)  were  made.  Com¬ 
plete  two  dimensional  data  were  taken  for  Ring  No.  1  but  only  the 
component  of  velocity  in  the  direction  of  ring  travel  was  mea¬ 
sured  for  Ring  No.  2. 

Two  samples  of  the  velocity  versus  time  data  from  Ring  No.  1 
are  shown  in  Fig.  22a  and  22b.  There  is  a  180°  ambiguity  in  the 
direction  of  velocity  (since  the  detected  frequency  shift  is 
always  positive)  which  is  easily  resolved  by  flow  visualization 
indicating  that  the  U  component  (Fig.  22a)  is  always  positive 
and  the  V  component  (Fig.  22b)  is  negative  to  the  left  of  the 
central  zero  and  positive  to  the  right.  Data  for  Ring  No.  2  is 
similar  with  velocities  being  five  times  greater  and  the  time 
scale  five  times  faster  since  the  ring  speed  UQ  of  Ring  No.  2 
is  about  five  times  that  of  Ring  No.  1. 
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Prior  to  data  analysis  the  Polaroid  oscilloscope  pictures 
are  photographically  enlarged  2.5  times  and  a  smooth  curve  drawn 
through  the  numerous  dots.  The  results  of  this  enlargement  are 
shown  in  Figs.  23  and  24.  These  curves  represent  a  portion  of 
the  data  taken  on  Ring  No.  1.  The  remaining  data  on  Ring  No.  1 
were  taken  in  the  vicinity  of  the  core  with  an  expanded  velocity 
scale. 

4.3.  Results  and  Discussion 

The  axial  component  of  velocity  U/UQ  versus  r/rQ  along 
z/rQ  =  0  for  the  two  vortex  rings  and  Hill's  spherical  vortex 
for  comparison  is  shown  in  Fig.  25.  The  value  of  rQ  is  defined 
by  the  zero  crossing  of  this  curve  and  the  core  radius  as  one-half 
the  distance  between  the  positive  and  negative  peak  velocities. 
This  gives  a  ratio  of  core  radius  to  ring  radius  of  a/rQ  =  0.27 
for  Ring  No.  1  and  a/rQ  =  0.075  for  Ring  No.  2,  indicating  that 
Ring  No.  2  is  a  relatively  thin  core  ring  while  Ring  No.  1  has  a 
relatively  thick  core  but  not  nearly  as  thick  as  Hill's  vortex 
where  a/rQ  -  1.414.  No  data  were  taken  in  the  core  of  Ring  No.  2 
due  to  the  lack  of  smoke  particles  in  the  vicinity  of  the  core. 

The  circulation  T(r)  is  obtained  by  integrating  the  velo¬ 
city  along  the  curve  C  (Fig.  21) ;  a  rectangle  with  one  side 
along  r  =  constant  and  the  other  sides  tending  to  infinity.  Thus 

r  =  (  «<*)- 

—  od 
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In  order  to  correct  for  the  finite  extent  of  the  data  it  is  assumed 

3 

that  (U (z) -U  ) ^1/z  after  the  last  data  point,  which  gives 

r  -  C'  («(*)-«.)  JZ  +  yzlu 

-«(?,■)  Z,]  ,5) 

where  z^  is  the  upstream  cutoff  and  is  the  downstream  cutoff. 

The  circulation  was  determined  by  finding  the  area  under 
curves  like  Fig.  22a  using  a  planimeter  and  adding  the  small 
correction  from  (5)  due  to  the  finite  extent  of  the  data.  This 
correction  was  always  less  than  5%.  The  total  circulation  is  the 
integral  along  r  =  0.  The  resulting  circulation  distribution 
(Fig.  26)  indicates  again  that  Ring  No.  2  is  the  thinner  core 
ring  having  the  circulation  more  concentrated  near  r/rQ  =  1.0. 

The  circulation  calculated  in  this  manner  can  give  some  indica¬ 
tion  of  the  vorticity  distribution  by  differentiating  the  curves 
of  Fig.  26.  However,  in  order  to  obtain  the  vorticity  this 
derivative  would  have  to  be  multiplied  by  an  appropriate  length 
corresponding  to  the  distribution  of  vorticity  in  the  z  direc¬ 
tion.  That  is,  the  slope  of  the  curves  of  Fig.  26  represents  the 
vorticity  weighted  by  an  axial  distribution  length  which  is 
unknown.  An  alternate  procedure  will  be  used  to  determine  the 
vorticity  distribution. 

A  summary  of  the  data  is  given  in  Table  2.  A  comparison 

between  the  actual  value  of  =  U  /T  /4ur  and  the  value  calcu- 

o  -  o  o 

lated  for  a  small  core,  solid  body  rotation  ring.  In  8r/a  -  1/4, 
shows  that  the  latter  is  30-40%  higher.  This  is  due  to  the  some¬ 
what  arbitrary  definition  of  a  and  the  fact  that  the  core  is 


not  in  solid  body  rotation  which  implies  A  has  a  value  other  than 


1/4. 


Since  both  components  of  velocity  were  measured  on  Ring 


No.  1  the  vorticity  distribution  throughout  the  vortex  ring  can 
be  found.  In  order  to  minimize  differentiation  of  data  the  fact 
that  to/r  is  constant  on  a  streamline  for  a  steady  axisymmetric 
flow  is  utilized.  The  streamlines  are  found  for  the  entire  flow 
field  and  the  value  of  u  is  calculated  only  along  z/rQ  =  0 
giving  a  value  of  w  on  each  streamline  so  w  is  known  every¬ 
where  . 

Integration  of  velocity  along  convenient  curves  gives  the 
stream  function  ip  so  that  the  streamlines ,  lines  of  constant 
can  be  plotted.  For  an  axisymmetric  flow 

d  ft  ~  rUo(r  -  rVdZ 

Integrating  from  z  =  +  along  a  line  r  =  constant  gives 

°°  r  ~z 

ft  (rscouvr.^  —  <LOViST,0°)  «  rj 


If  we  take  V  ^  1/z  after  some  cutoff  z^  and  non-dimensionali- 


zing 


4.  -L  =  jl  v 

r0l  Uo  *  vx./  To  i-  t.U0r0 

-  c 


The  stream  function  is  thus  determined  by  integrating  curves 
like  Fig.  22b,  with  the  results  shown  in  Fig.  27. 
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As  a  check  the  stream  function  along  z/rQ  =  0  is  calculated 
by  integrating  the  U  component  of  velocity  in  Fig.  25  using 


^(r,o) 
r„l  u0 


u  ( r/r.l 


u, 


The  values  of  ^  along  z/rQ  =  0  from  the  two  methods  are 
shown  in  Fig.  28.  This  critical  check  of  calculating  ip  from 
different  data  sets  indicates  that  the  overall  consistency  of  the 
data  is  very  good.  The  actual  streamlines  are  shown  superposed 
on  a  strobe  photograph  of  Ring  No.  1  in  Fig.  29. 

The  vorticity,  w,  is  calculated  along  z  =  0  differentiating 
the  data 


cuCr>  o) 

(«»/0 


Hu1.) 


«-o 


The  first  term  is  the  slope  of  the  curves  of  Fig.  22b  at  z  =  0 
and  the  second  term  is  obtained  by  differentiating  the  smoothed 
curve  of  Fig.  25.  The  results  (shown  in  Fig.  30)  of  differen¬ 
tiating  the  first  term  show  scatter  in  the  immediate  vicinity  of 
r/rQ  =  1  caused  in  part  by  the  directional  ambiguity  of  the  LDV 
system.  The  V  component  of  velocity  rapidly  changes  from  a  large 
positive  value  to  a  large  negative  value  giving  data  curves 
(Fig.  22b)  in  which  the  determination  of  the  slope  at  z  =  0  is 
difficult.  The  vorticity  distribution  along  z  =  0  (Fig.  31) 
shows  the  expected  concentration  of  vorticity  in  the  core. 

Using  the  fact  that  ui/r  is  constant  on  streamlines  in  con¬ 
junction  with  Fig.  29  and  31,  we  can  find  the  vorticity  at  any 
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point  in  the  ring.  it  is  estimated  that  roughly  85%  of  the  vor- 

2 

ticity  is  within  the  ^/rQ  UQ  -  • 4  streamline. 

From  the  circulation  distribution  (Fig.  26)  it  is  expected 
that  Ring  No.  2  would  have  an  even  more  concentrated  vorticity 
distribution.  These  results  are  in  contrast  to  Maxworthy's  (Ref  30) 
claims  that  the  vorticity  is  more  uniformly  distributed  through¬ 
out  the  vortex  ring.  It  is  felt  that  this  discrepancy  is  mainly 
due  to  the  inaccuracy  of  the  hydrogen  bubble  technique  although 
Maxworthy's  rings  may  be  slightly  thicker  than  the  ones  pre¬ 
sented  here. 

Also  shown  in  Fig.  31  are  the  vorticity  distribution  for  a 
small  core  solid  body  rotation  ring  (using  the  measured  value  of 
a/rQ)  and  Hill's  spherical  vortex.  This  comparison  shows  that 
for  rings  of  the  same  size  and  velocity  of  translation  neither 
model  properly  describes  the  experimentally  produced  vortex  rings. 

In  a  steady  axisymmetric  flow  u>/r  =  F ( ip)  .  Most  theories 
assume  F(i|>)is  a  constant  since  this  is  the  only  tractable  case, 
but  Fig.  32  shows  that  this  is  not  the  case. 

4.4.  Ring  Properties  from  Vorticity  Measurements 

From  the  measured  vorticity  distribution  along  z/rQ  =  0, 
we  can  calculate  the  circulation  and  the  velocity  parameter  ^ 
by  using  the  theory  presented  by  Bliss  (Ref  24).  In  this  work, 
the  solution  for  the  vorticity  distribution  within  a  slightly 
curved  vortex  filament  was  obtained  by  the  method  of  matched 
asymptotic  expansions.  This  theory  shows  that  the  vorticity 
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distribution  in  a  slightly  curved  filament  may  be  written 

to(r)-  f  SIU© 

where  uq  is  the  vorticity  distributed  in  a  straight  filament  with 
radial  symmetry  and  is  the  small  perturbation  due  to  curvature, 
uj  is  symmetric  about  r  =  rQ  (P  =  0)  and  is  antisymmetric. 

We  may  use  these  properties  to  obtain  u)q  directly  from  the  mea¬ 
sured  data  since  w  is  then  the  average  value  of  the  measured 
vorticity  for  points  equally  spaced  about  r/rQ  =  1. 

If  uo  is  known  for  a  vortex  ring,  the  circulation  and  ring 

velocity  can  be  calculated.  The  functions  w  and  a).  are  shown 

o  1 

in  Fig.  23. 

The  circulation  of  the  ring  in  terms  of  w  is  just 

Direct  numerical  integration  of  u)Q  calculated  from  the  measured 

2 

vorticity  w  gives  a  value  for  r0  of  1.17  ft  /sec  in  close 

2 

agreement  with  the  measured  value  1.21  ft  /sec  (Table  1). 

The  non-dimensional  ring  velocity  ^  can  also  be  determined 
from  the  measured  vorticity.  The  non-dimensional  self-induced 
velocity  of  a  vortex  ring  with  an  arbitrary  distribution  of 
vorticity  is  given  by 

v  =  Jtn  SJGl  -4-  K-Vi.  =  8  -  +  ( A “ ^ a/rb) 
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where 


A  =  Mm  f  -~?np 

p  ao  '-to  ' 


v  is  the  non-dimensional  swirl  velocity  (in  p,0  coordinate  sys- 

A 

tem)  v  =  v  2ra/r  and  p  =  where  a  is  some  typical  dimen- 
o  o  o  r 

o 

sion  in  the  vortex  core.  (Actually,  the  parameter  A  -  In  a/rQ 
is  independent  of  the  choice  of  length  scale.) 

Straightforward  manipulation  yields  a  formula  for  the  group 
A  -  In  a/ro  in  terms  of  the  zeroth-order  vorticity  ujo.  The 


results  being  _  r  -  W.(^)^n/5  £  f'  C00{f')JLp'd(.f 

A  Vy,  -  -is - - - ° - - - - 
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If  the  vorticity  is  uniform  within  a  radius  r  =  a,  this  expres¬ 
sion  gives  the  familiar  result  A  =  1/4.  Direct  numerical 
integration  of  this  formula  using  the  values  of  u jq  obtained 
from  the  measured  vorticity  gives  ^  =  3.01  as  compared  with  a 
measured  value  of  2.46,  an  error  of  20%.  This  is  in  somewhat  better 
agreement  than  the  value  3.14  obtained  by  assuming  uniform  vor¬ 
ticity.  Since  the  ratio  of  core  size  to  radius  is  0.27  perhaps 
a  20%  discrepancy  between  experiment  and  an  asymptotic  theory 
which  assumes  a/r  <<  1  is  not  surprising. 

A  value  of  A  =  0.136  is  obtained  for  the  vorticity  distribu¬ 
tion  of  the  ring  for  which  a/ rQ  =  0.27. 

4.5.  Stability 

During  the  study  of  the  structure  of  vortex  rings  it  was 
observed  that  rings  were  unstable  to  azimuthal  perturbations,  a 
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fact  that  had  been  reported  by  Krutzsch  (Ref  18)  in  1939.  An 
analytical  study  by  Widnall  (Ref  42)  concluded  that  all  vortex 
rings  were  unstable  to  these  azimuthal  perturbations  with  the 
number  of  unstable  waves,  n,  being  a  function  of  the  non- 
dimensional  ring  velocity  ^  =  \JQ/r Q/4vrQ.  For  a  given  vortex 
core  only  one  mode  is  unstable  and  the  smaller  the  vortex  core 
(i.e.  the  higher  the  value  of  ^)  the  larger  the  number  of  waves 
in  the  mode.  This  result  is  shown  in  Fig.  33  from  Ref  42. 

A  side  view  multiple-flash  strobe  photograph  of  an  n  =  7 
mode  instability,  Fig.  34,  shows  that  the  ring  appears  to  be 
completely  destroyed.  For  a  higher  mode  case,  another  ring  travel¬ 
ing  at  a  slower  velocity  is  formed  after  the  initial  instability. 

Front  view  photographs  of  an  n  *  6  and  n  =  7  mode  insta¬ 
bility  are  shown  in  Figs.  35  and  36.  The  actual  observed  mode 
number  n  of  the  instability  is  shown  in  Fig.  37  as  a  function 

of  =  U  /r  /4Tir  .  The  LDV  system  was  used  to  determine  V  ,  and 
o  o  o  o 

UQ  and  rQ  were  measured  photographically.  Also  shown  is  the 
value  of  n  that  is  predicted  to  be  the  most  unstable  for  that 
value  of  V.  Although  the  trend  of  increasing  n  with  increasing 
(decreasing  core  size)  is  in  agreement  with  the  theoretical 
predictions,  the  actual  values  of  n  vs  ^  are  not  in  good 
agreement.  This  is  most  likely  due  to  the  fact  that  the  ratio 
of  wave  length  to  vortex  core  radius  of  the  unstable  mode  is  not 
small  as  required  by  the  asymptotic  theory.  The  measured  radial 
perturbation  of  the  ring  due  to  the  instability  as  a  function  of 
time  is  shown  in  Fig.  38.  The  growth  of  the  perturbation  during 
the  early  stages  of  the  instability  shows  good  agreement  with  the 
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CHAPTER  V 

HELICOPTER  ROTOR  WAKE 

5.1.  Introduction 

The  wake  of  a  helicopter  rotor  is  another  example  of  a  flow 
field  where  the  use  of  conventional  instrumentation  is  hampered 
by  probe  interference  and  calibration  in  a  rapidly  changing 
velocity  field.  Hollow  rotor  blades  were  constructed  so  that  smoke 
could  be  injected  at  any  point  along  the  blade  to  provide  flow 
visualization  and  particles  for  the  LDV.  From  the  velocity  iea- 
surements  the  bound  circulation  and  hence  the  load  distribution 
along  the  rotor  are  found. 

5.2.  Construction  of  a  Rotor  for  Use  with  LDV  System 

A  rotor  stand  capable  of  spinning  a  one  or  two  bladed  model 
motor  was  designed  and  constructed.  The  hollow  rotor  blades  per¬ 
mit  smoke  to  be  deposited  at  the  tip  or  any  other  desirable 
place  along  the  blade. 

A  cross  section  of  the  blade  in  Fig.  39  shows  the  epoxy¬ 
fiberglass  structure  with  a  stainless  steel  spar,  beryllium- 
copper  inner  chamber  and  a  tungsten  rod  to  maintain  balance  about 
the  quarter  chord.  Initially  an  aluminum  model  was  machined  from 
which  a  mold  (Fig.  40)  was  cast  of  Stycast  2850,  a  high  strength 
casting  resin. 

The  beryllium-copper  was  spot-welded  to  the  stainless  steel 
spar  and  the  tungsten  rod  silver  soldered  on  the  front  of  it.  An 
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inner  mandrel  was  inserted  inside  the  beryllium- copper  and  the 
outside  wrapped  with  fiberglass  wet  with  Eccomold  L-28  lamina¬ 
ting  resin.  This  structure  was  compressed  in  the  mold  to  form 
the  NACA-0012  airfoil  section.  After  the  epoxy  cured  the  blade 
was  removed  from  the  mold  and  the  inner  mandrel  removed.  The 
resulting  hollow  blades  are  shown  in  Fig.  40b  and  40c. 

The  drive  stand  for  the  rotor  (Fig.  41)  has  a  hollow  shaft 
to  allow  smoke  to  pass  to  the  rotor.  The  shaft  is  driven  by  an 
1800  rpm  synchronous  motor. 

Table  3  lists  the  relevant  parameters  for  the  single  bladed 
rotor  used  throughout  this  experiment. 

5.3.  Flow  Visualization  of  Rotor  Wake 

Several  interesting  features  of  the  flow  field  in  the  wake 
of  the  model  rotor  were  observed  using  smoke  flow  visualization. 
The  coordinate  system  used  for  describing  the  rotor  experiments 
is  shown  in  Fig.  42. 

A  side  view  strobe  picture  of  the  tip  vortex  (Fig.  43) 
shows  an  unexplained  wavelike  structure.  Oscillations  in  the 
tip  vortex  of  a  single  bladed  rotor  are  mentioned  by  Ref.  43 
but  no  photographs  or  explanation  are  given.  Attempts  have  been 
made  to  explain  this  wave  using  the  analysis  of  the  stability  of 
a  helical  vortex  by  Widnall  (Ref.  44) ,  but  the  phase  of  the 
observed  wave  is  different  than  the  theory  predicts.  This  is 
probably  due  to  the  fact  that  the  theory  does  not  include  the 
effects  of  the  inboard  vortex  sheet. 

When  smoke  is  released  upstream  of  the  rotor  plane,  it 
becomes  entrained  in  the  core  of  the  tip  vortex  and  travels  up 
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the  core  towards  the  rotor  blade  as  shown  in  Fig.  4  4  indicating 
a  strong  axial  flow  in  the  vortex.  Since  flow  along  the  core  of 
a  vortex  affects  the  self-induced  motion  this  may  be  important  in 
analytic  determination  of  vortex  position.  A  slit  was  cut  in  the 
top  of  the  blade  at  the  quarter  chord  between  r/Ro  =  . 82  and 

r/Ro  =  .94.  The  results  of  releasing  smoke  through  this  slit  and 

at  the  tip  (Fig.  45) ,  show  that  the  sheet  of  smoke  does  not 

immediately  roll  up  into  the  tip  vortex.  Additional  observations 

indicate  that  it  takes  over  a  complete  revolution  before  the 
sheet  is  rolled  up. 

In  order  to  obtain  the  position  of  the  tip  vortex,  smoke 
was  ejected  at  the  tip  of  the  blade  and  the  wake  illuminated  with 
a  narrow  plane  of  light.  From  a  series  of  photographs,  examples 
of  which  are  shown  in  Fig.  46,  the  tip  vortex  position  is  estab¬ 
lished.  The  coordinates  of  the  tip  vortex  given  in  Figs.  47  and 
48  are  in  general  agreement  with  the  averaged  experimental  re¬ 
sults  of  Ref.  43. 

A  photograph  looking  along  the  positive  z  axis  from  beneath 
the  rotor,  Fig.  49,  shows  some  interesting  features  of  the  tip 
vortex.  It  can  be  seen  that  the  tip  vortex  gradually  increases 
in  diameter  until  regions  of  sudden  enlargement  (vortex  break¬ 
down)  are  encountered.  A  graph  of  core  diameter  vs.  azimuthal 
angle.  Fig.  50,  shows  this  increase. 

The  two  component  laser  Doppler  velocimeter  (LDV)  des¬ 
cribed  in  Chapter  II  v/as  used  to  measure  velocities  around  the 
hovering  model  rotor  blade.  Two  separate  measurement  setups 
were  used  to  determine  all  three  components  of  the  velocity 
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vector.  Several  examples  of  the  velocity  data  are  shown  in 
Fig.  51  for  a  single  untwisted  blade  at  an  anc>le  of  attack 
0O  =  10°.  The  curves  show  velocity  vs.  time  as  the  rotor  blade 
passes  a  fixed  point.  The  sign  of  the  velocity  component  is 
not  given  by  the  LDV  system,  due  to  an  inherent  180°  ambiguity, 
and  must  be  determined  by  independent  means  such  as  flow  visuali¬ 
zation.  The  proper  sign  of  the  velocity  component  is  indicated 
for  each  data  curve.  Fig.  51  a-f  are  on  a  short  enough  time 
scale  so  that  details  of  the  flow  field  in  the  immediate  vicinity 

of  the  blade  are  shown.  Fig.  51  e  shows  the  downwash  v  for  a 

z 

complete  revolution  of  the  blade.  Note  that  the  data  below  the 
blade  (Figs  51b, d,  and  f)  in  the  relatively  small  region  of  the 
viscous  wake  are  extremely  unrepeatable. 

By  integrating  the  velocity  data  along  the  contour  shown 
in  the  dotted  lines  in  Fig.  42,  the  bound  circulation  and  hence 
the  section  load  can  be  determined.  The  entire  load  distribu¬ 
tion  is  found  by  performing  the  integration  at  a  series  of  radial 
locations.  The  results  of  these  integrations  are  shown  in 
Figs.  5 2  and  53.  From  blade  element  theory  (Ref.  45)  for  an 
untwisted  rotor 
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and  J[  C-f 

f  PM\ 

fix 

*  VfMOTTJ V 

where 

T  (x)  =  bound  circulation  at  x 
x  =  r/Ro 

Ro  =  radius  of  rotor 
fl  *  angular  velocity 
o  solidity 
0#  “  angle  of  attack 
a  =  arve 

CT  =  thrust  nt 

The  comparison  between  blade  element  theory  and  experiment  shows 
the  expected  discrepancy  at  the  tip  since  the  theory  does  not 
include  the  effect  of  the  tip  vortex.  When  a  tip  loss  factor 
of  .96  is  applied  to  the  theory,  a  value  of  CT  =  .00388  is 
obtained  which  is  within  1.8%  of  the  value  obtained  from  an 
overall  force  balance.  A  value  of  CT  =*  .00357  is  obtained  from 
the  laser  data  by  finding  the  area  under  the  experimental  curve 
of  Fig.  53.  Although  this  is  in  error  by  6.5%,  increased  accu¬ 
racy  could  be  obtained  by  taking  more  data  and  averaging  and  by 
improving  the  accuracy  of  the  Doppler  frequency  measurement. 

Also  the  method  of  obtaining  the  local  lift  from  pU(x)T(x)  is 
in  error  in  the  tip  region.  The  load  distribution  of  Fig.  53 
and  the  tip  vortex  position  of  Figs.  47  and  48  offers  a  relatively 
simple  test  case  for  any  free  vortex  wake  analysis  of  a  hovering 
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5.4.  LDV  Measurements  Near  the  Tip  Vortex 

Numerous  attempts  were  made  to  obtain  data  in  the  tip 
vortex  of  the  model  rotor.  Examples  of  LDV  data  taken  at  a 
position  strongly  influenced  by  but  not  directly  within  the 
tip  vortex  are  shown  in  Fig.  54.  The  point  of  closest  approach 
of  the  vortex  to  the  measuring  volume  occurs  at  approximately 
50°  behind  the  quarter  chord.  The  data  are  now  beginning  to  show 
more  scatter.  When  the  probe  is  moved  closer  to  the  tip  vortex 
the  scatter  renders  the  data  uselsss.  Since  the  LDV  data  near  the 
blade  are  repeatable  the  cause  of  scatter  must  be  due  to  fluctua¬ 
tions  of  the  tip  vortex  position.  In  order  to  obtain  information 
about  the  structure  of  the  tip  vortex  an  electronics  system  that 
does  not  require  the  flow  field  to  be  periodic  is  needed.  It  is 
also  recommended  that  a  larger  model  rotor  be  used  so  that  the 
reynolds  number  would  be  higher  and  the  velocity  gradients  across 
the  probe  smaller. 
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CHAPTER  VI 

CONCLUSIONS  AND  RECOMMENDATIONS 

6.1.  Summary  and  Conclusions 

The  use  of  conventional  instrumentation  such  as  a  hot  wire 
anemometer  for  measuring  velocity  in  vortex  filament  is  hampered 
by  probe  interference  and  the  difficulty  of  calibrating  the  probe 
in  a  flow  field  where  the  velocity  is  changing  rapidly  in  magni¬ 
tude  and  direction.  For  these  reasons  a  laser  Doppler  veloci- 
meter  (LDV)  is  used  throughout  these  experiments.  The  LDV 
measures  the  velocity  of  smoke  particles  by  measuring  the 
Doppler  shift  of  laser  light  scattered  from  the  particles.  The 
LDV  system  designed  and  constructed  is  a  dual  scatter  system 
capable  of  determining  two  components  of  the  velocity  vector. 

The  polarization  of  the  incident  beams  used  for  one  component 
is  perpendicular  to  that  for  the  other  component  to  prevent  cross 
talk  and  allow  the  components  to  be  separated.  By  using  conven¬ 
tional  prism-type  beamsplitters  a  compact  in  line  system  that 
maintains  equal  path  lengths  for  all  beams  has  been  constructed. 

The  most  unique  feature  of  the  LDV  system  is  the  electronics  sys¬ 
tem  which  is  capable  of  handling  the  large,  fast  changes  in  velocity 
which  are  encountered  in  the  flow  fields  of  interest.  The  electron¬ 
ics  system  described  in  section  2.5  requires  that  the  flow  field  be 
periodic.  It  is  a  "sampling"  system  in  that  only  a  few  points  at  a 
particular  phase  on  the  velocity-versus-time  curve  are  determined 
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per  cycle.  By  varying  the  phase  the  remaining  points  are  deter¬ 
mined.  The  system  requires  only  a  spectrum  analyzer  for  demodu¬ 
lation  and  a  storage  oscilloscope  for  retaining  the  data. 

The  main  drawback  to  the  LDV  is  the  fact  that  the  velocity 
measured  is  that  of  smoke  particles  and  not  the  fluid  velocity. 
Therefore  it  is  necessary  to  assure  that  the  particle  size  is 
small  enough  so  its  velocity  is  nearly  equal  to  the  fluid  velo¬ 
city.  This  is  especially  true  in  the  present  experiments  where 

3 

oil  particles  with  a  density  of  10  times  that  of  air  are  used. 
The  oil  particles  are  generated  by  a  standard  mist  lubricator — a 
device  normally  used  to  lubricate  high  speed  bearings.  The 
particles  have  a  diameter  of  approximately  1.0/x  which  was  deter¬ 
mined  by  a  light  scattering  experiment  that  measured  the  shape 
of  the  forward  scattering  lobe.  It  is  estimated  that  oil  par¬ 
ticles  of  this  size  will  follow  the  fluid  velocities  encountered 
in  these  experiments  within  1.0%. 

A  detailed  study  of  the  structure  and  stability  of  vortex 
rings  using  the  two- component  LDV  system  has  been  made.  Vortex 
rings  were  generated  by  pulsing  air  through  a  sharp-edged 
orifice  using  a  loud  speaker.  The  circulation  and  core  size  are 
controlled  by  the  amplitude  and  duration  of  the  electric  signal 
applied  to  the  loud  speaker. 

Detailed  surveys  of  two  vortex  rings  were  made;  a  rela¬ 
tively  thick  core  ring  with  a  ratio  of  core  radius  to  ring 
radius  a/rQ  =  *27 ,  and  a  relatively  thin  core  ring  witji 
a/rQ  =  .075.  The  principal  result,  Fig. 31, shows  that  the  vorti- 
city  is  concentrated  in  a  small  core  region  but  that  the  core  is 
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not  in  solid  body  rotation.  The  self-induced  velocity  predicted 
by  theory  for  the  measured  vorticity  distribution  is  in  error  by 
20%,  probably  due  to  the  fact  that  a/rQ  3  .27  is  not  small  enough 
for  the  theory  to  be  valid.  (Only  the  thick  core  case  could  be 
compared. ) 

During  the  course  of  the  investigation,  the  fact  that  vor¬ 
tex  rings  are  unstable  was  rediscovered.  Azimuthal  perturba¬ 
tions  with  between  6  and  16  waves  spontaneously  develop  on  the 
perimeter  of  the  ring  and  continue  to  grow.  For  a  given  vortex 
core,  only  one  mode  is  unstable  and  the  smaller  the  vortex  core, 
the  larger  the  number  of  waves  in  the  mode.  Theoretical  pre¬ 
dictions  of  the  mode  number  are  not  good,  but  the  correct  trend 
is  indicated.  The  amplification  rate  is  accurately  predicted 
probably  because  all  modes  have  about  the  same  amplification  rate. 

A  single  bladed  model  helicopter  rotor  has  been  constructed. 
The  blade  is  hollow  to  allow  smoke  to  be  deposited  at  the  tip  or 
any  desirable  point  along  the  span  providing  flow  visualization 
and  scattering  particles  for  the  LDV  system.  The  flow  visuali¬ 
zation  data  in  hover  show  an  unexplained  wave like  structure  on 
the  helical  tip  vortex.  Also  a  significant  axial  velocity  up 
the  core  of  the  tip  vortex  towards  the  blade  is  evident  from 
the  smoke  photographs. 

Velocities  around  the  hovering  blade  were  measured  with 
the  LDV  system.  By  making  two  separate  measurement  setups,  all 
three  components  of  the  velocity  vector  were  determined.  .  From 
this  data  the  bound  circulation  and  hence  the  spanwise  load 
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distribution  was  found.  Attempts  at  using  the  LDV  system  to 
measure  the  structure  of  the  tip  vortex  were  largely  unsuccessful 
due  to  random  fluctuations  in  the  position  of  the  vortex  which 
could  not  be  handled  with  the  present  electronics  system. 

6.2.  Recommendations 

In  order  to  improve  velocity  measurements  in  vortices  within 
an  LDV  it  is  first  necessary  to  develop  fast  electronics  that  can 
handle  the  large  and  rapid  changes  in  Doppler  frequency  that  are 
encountered.  Although  the  presently  available  frequency  trackers 
are  not  fast  enough,  their  response  to  large  fluctuations  could 
be  improved  by  increasing  the  I.F.  frequency.  If  the  I.F.  fre¬ 
quency  is  made  much  larger  than  the  changes  in  Doppler  frequency 
then  the  time  response  should  be  improved. 

Once  adequate  electronics  is  developed  it  is  necessary  to 
overcome  the  problem  of  random  fluctuations  in  the  position  of 
the  vortex.  This  could  be  done  by  rapidly  moving  the  measuring 
volume  of  the  LDV  through  the  vortex.  This  scan  must  be  fast 
compared  to  the  fluctuation  speed  of  the  vortex  so  that  the  vor¬ 
tex  position  is  essentially  steady  during  the  measurement.  Ulti¬ 
mately  what  is  needed  is  a  three-dimensional  raster  scan  of  the 
volume  containing  the  vortex.  While  the  method  of  producing 
such  a  scan  is  not  completely  clear,  it  would  appear  that  com¬ 
mercially  available  pizeoelectric  or  electromagnetic  scanners 
should  be  adequate. 

The  study  of  vortex  rings  could  be  expanded  using  the  LDV 
setup  without  modification.  An  investigation  of  the  rollup  and 
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formation  process  would  be  particularly  interesting.  The  pre¬ 
sent  investigation  could  also  be  extended  to  include  rings  with 
smaller  cores, and  the  effects  of  buoyancy. 

The  use  of  the  improved  LDV  system  described  above  for 
measurements  of  the  flow  field  around  the  helicopter  rotor  is 
the  most  practical  and  pressing  problem.  The  use  of  a  larger 
model  rotor  would  increase  the  Reynolds  number  and  also  the  size 
of  the  tip  vortex  compared  to  the  probe  volume.  With  this  sys¬ 
tem  the  structure  and  position  of  the  tip  vortex  and  the  load 
distribution  could  be  found. 
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TABLE  1 


LDV  Characteristics 

Laser  -  HeNe  Spectra  Physics  Model  124 

A  =  632  8A 

P  =  15  raw 
Li 

2b  =  3.8  mm  (at  focusing  lens) 
f  =  12  in.  =  304.8mm  (focusing  lens) 

0.72  =  .0681  Rad.  =  3.90° 

3 

fD  =  (6.5716  x  lO^U)Hz  (where  U  is  the  velocity  in  ft/sec) 

b  -  =  .032  mm  =  .0013" 

o  *TT 

Diameter  of  probe  volume  =  2bQ  =  .0026" 

Length  of  probe  volume  =  2bQ/sin  0fi/2  =  .038" 

-6  3 

Volume  of  probe  =  2.02  x  10  cm 
Particles  are  made  from  "Fog  Juice"  (light  oil) 
a  =  Tid/A  =  5.0 
N  =  8.5  x  10^  particles/cm"* 


TABLE  2  DATA  SUMMARY 
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Table  3 

Model  Rotor  Parameters  (Single  Blade) 

C  Chord  =  1.00" 

R  =  7.203" 

o 

9  =  10° 

o 

fl  =  1800  RPM 

ARq  =  tip  speed  *  113.1  ft. /sec. 

<T  =  .0442 

flR  C 

Re  =  =  62,830 

CT  *  .00381  (Measured  by  Force  Balance) 

CT.  =  .0862 

' a 


FREQUENCY  (Hi) 


VELOCITY  (ft/sec) 

Fig.  2.  Doppler  shift  vs.  velocity  for  He-Ne  laser  at  various 
scattering  angles. 
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REFERENCE  BEAM  SYSTEM 
Fig.  3a 


DUAL  SCATTER  SYSTEM 
Fig.  3b  * 


Fig.  4a.  Doppler 
signal  from  .0002 
dia.  wire  passing 
through  center  of 
probe  volume. 

Vert. scale  lOOMV/cm 
Horiz.  "  . lMS/cm 


Fig.  4b.  Doppler 
signal  from  .0002 
dia.  wire  which  did 
not  pass  through 
center  of  probe  vol. 
Scale  same  as  4a. 


Fig.  4c.  Doppler 
signal  from  spin¬ 
ning  disk. 

Vert. scale  lOOMV/cm 
Horiz.  "  5C*a/cm 
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Dual-Scatter  System 


Demodulation  of  Test 

£  -  500  kHz 

o 

Df  -  450  kHz 
f  -  33  cy/sec 


ZHH)  0 


Vortex  ring. 


U  VELOCITY 
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Velocity  distribution  through  vortex  ring. 
Time  Scale :  5  nis/cm. 
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Fig.  20.  Scheme  of  Particle  Sizing  Based  on 
Scattering  Measurements  Within  the 
Forward  Lobe.  The  Solid  Curves  are 


Ratios  of  l,  versus  a.  i.e.  the  10° 
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Vortex  Ring 


Fig.  21.  Coordinate  system  for  vortex  ring  data  analysis. 


Fig.  22a.  U  component  of  velocity 
at  r/rQ  =-. 876  for 

Ring  No.  1. 

Time  Scale  20  MS/cm 
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Fig.  22  b.  V  component  of  velocity 
at  r/r  =  -  1.08  for 
Ring  No.  1. 

Time  Scale  20  MS/cm 
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Fig.  26.  r/r  versus  r/r 
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Fig.  28.  */rQ2U0  versus  r/rQ  at  Z/rQ  *  0. 


^  -  Calculated  fron  Eq.  2. 

Solid  curvea  calculated  fron  Eq.  3. 


Fig.  29. 


Side  view  of  ring  No.  1  with  measured 
stream  lines  superposed. 


Fig.  30.  Partial  Derivatives  Used  For 
Calculating  the  Vortioity. 
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Fig.  37.  Number  of  waves  in  the  unstable  mode  as  a  function  of 
X  experiment?  0  theoretical  prediction  at  maximum 
amplification. 


Fig.  39-  Cross  Section  of  Blade  -  NACA  0012. 
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Fig.  43.  Tip  vortex  showing  wave  structure 


Fig.  44.  Strobe  photograph  showing  flow  up  the  tip 
vortex  towards  the  blade. 
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Position  of  Tip  Vortex 


Radial  and  axial  position  of  tip  vortex. 


Strobe  photograph  showing  a  gradual  increase 
in  core  diameter  and  vortex  breakdown. 
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Fig.  50.  Cere  diameter  versus  azimuthal  angle 
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Bound  Circulation  versus  Radius 


Nondimensional  Load  Distribution  Along  Blade 


